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Abstract. Jet reconstruction in heavy-ion collisions is strongly affected by soft 
background from the underlying event. For an appropriate interpretation of the 
jet observables it is essential to understand the influence of the background and its 
fluctuations on the reconstructed jets. With this purpose we study random cones 
and the response of a known probe embedded in a heavy-ion events. The embedded 
probe can be a single high-px track or a jet from a simulated or real pp event. This 
allows a detailed study of background fluctuations and veriflcation of the performance 
of background subtraction methods. 



1. Introduction 

We present our results on region-to-region energy fluctuations of the soft background in 
Pb-Pb collisions at the LHC (a/snn = 2.76 TeV) measured with the ALICE experiment. 
The studies are based on charged tracks with a low px cut of 150MeV/c. We used two 
different methods random cones and fast embedding to investigate the background. For 
fast embedding we discuss the case of embedded single high-pT tracks. 

A good knowledge of the background fluctuations is mandatory for the 
interpretation of jet measurements. Recently, for example the impact of background 
fluctuations is of particular interest since the flrst measurements of imbalanced jets in 
Pb-Pb collisions at the LHC have been published [l] [2] [s]. 

2. Background Subtraction 

The average background has to be estimated and subtracted on event-by-event basis. 
For this purpose, the heavy-ion event is clusterized with the kx-algorithm from the 
Fast Jet package 14]. We assume that most of these clustered objects consist of soft 
background, even though this is not the case for all of them. We call these background 
jets in the following. The clustering procedure allows to eliminate regions with strong 
deviations from the average background, as we determine the background density p of 
the event as the median of those background jets: 



p = median(pJp/A!.gJ, 
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where is the momentum of the background jet i with the jet area A^g^ in the rj — (j) 
plane. For the median we exclude the two hardest jets. In figure [T]p is shown as function 
of centrality (a) and multiplicity (b). The background density increases linearly with 
the raw number of track^ (multiplicity). 




Figure 1: Background density p in Pb-Pb collisions at ^/s = 2.76 TeV as function of 
centrality (a) and multiplicity (b) with low track cut off at 150MeV/c. Multiplicity 
is in this case the raw (uncorrected) number of tracks in the event, which is the number 
of input tracks for the jet finder. 

p X A jet is the expected amount of background energy which contributes to the jet 
momentum. For a jet radius of i? = 0.4 and corresponding jet area of Aj^t = 0.5 
it is above 50GeV/c and up to lOOGeV/c for the 10% most central events. The 
average background p is estimated event-by-event and p x Aj^t is subtracted from the 
reconstructed jet momentum. What still persists are the region-to-region background 
fiuctuations 6pt:, which can not be corrected event-by-event and needs to be convoluted. 
The jet energy resolution is limited by these fiuctuations js]: 

PT,jet = PTjet " P X A jet " ^PT- 

3. Jet Reconstruction 

The /ct algorithm [6] [7j sequentially combines the tracks i and j with the smallest 
distance parameter dij = min(p^ j,p^j)— where Ai?jj = J (Arj)'^ + (A^)^ and R is 
the radius parameter. The algorithm is especially convenient for background jets. 
Since the dij is weighted with the squared transverse momentum PT(ij)' basically 

I Raw number of reconstructed charged tracks after quality cuts, which is the number of input tracks 
for the jet finder. 
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starts clustering with low-px tracks and it accordingly is adaptable to soft contributions. 
This is disadvantageous for the reconstruction of real jets in an environment of large 
soft background (see e.g. [s], (sj). 

The anti-^T algorithm on the contrary weights the distance between the tracks 
with the inverse transverse momentum min(l/p^ •, l/p\ it accordingly mainly starts 
clustering with high-pT particles and is much more robust against soft background in 
terms of area stability and back reaction [s]. Hence, we use the anti-fcx algorithm for 
real jet reconstruction in Pb-Pb collisions. The jets are reconstructed with R = 0.4 in 
a pseudorapidity range according to the jet axis of |?7| < 0.4. 

4. Fast Embedding 

One of the methods we use to determine the background fluctuations is the embedding 
of a known probe into a real heavy- ion event (see e.g. |9j). The probe can be single 
high-pT tracks or fully simulated (e.g. PYTHIA + GEANT) or real pp jet events. Here 
we discuss only the results from embedded single-tracks, which are taken as jets with 
only one (high-px) track. This method is called fast embedding, since the embedding 
is done for already reconstructed tracks and does not take into account possible effects 
on reconstruction, like possible merging of track clusters with influence on the tracking 
efficiency and resolution. 

Since we know the transverse momentum of the embedded probe we can calculate 
the background fluctuations 6pT after background subtraction and subtraction of the 
momentum of the embedded probe [o] : 

That allows us to study in detail the influence of the soft background on the 
reconstructed jet observables. We also can verify the performance of the background 
subtraction methods we use. If the average background is correctly subtracted the mean 
of the distribution is around zero, otherwise a systematic shift in Sp^ will occur. 

A certain matching of the embedded probe and the reconstructed jet in the heavy- 
ion event is necessary. For embedded single tracks we just match the track with the 
reconstructed jet which contains the track. As additional requirement both jets (probe 
and reconstructed) need to be in the jet acceptance. 

The single tracks are randomly embedded over full jet r/ — acceptance (|?7| < 0.4) 
with a flat pt distribution from 50GeV/c to 250GeV/c. We need to keep in mind 
that the probes are embedded in a flat centrality distribution, that does not equate to a 
corresponding centrality distribution of an inclusive jet spectrum, since the cross-section 
of jet production in a heavy-ion event increases with the number of binary collisions. 

5. Random Cones 

An alternative method to study the background of heavy-ion events are random cones. 
Here cones with an fixed size of i? = 0.4 are randomly placed in the heavy- ion event 
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and the in the cones are summed up to background jets. The average background is 
subtracted as described earher, so we directly get the background fluctuations: 



6. Results 

Figure [2] shows the background fluctuations 5pT for two centrality classes comparison the 
two discussed methods random cones and fast embedding of single tracks. Figure 
shows the results for the 10 % most central events. Since the Sp^ distribution is 
mainly caused by uncorrected Poissonian fluctuations it is supposed to be a gamma- 



distribution 10 . However, for a large number of input tracks, which is the case for the 
most central events, the distribution can be well described with a Gauss in ideal case of 
only uncorrelated fluctuations. But, in addition jets are already part of the heavy-ion 
events, and those cause the tail on the right-hand side of the 6pT distribution. So we 
specify the distribution with a Gaussian fit on the left-hand side as a first assumption. 
More precisely it is a iterative fit with a fit range from fi — 3a to /i + 0.5cr with the 
mean fi and width a of the Gaussian fit. The termination condition is a shift of /i less 
than 0.1 GeV/c, the maximum number of iterations is 20. 

The results from random cones and track embedding agree pretty well. For both 
the mean is around zero, that demonstrates that the background subtraction method 
works as desired. The results demonstrate that the impact of the average background 
in the event to the reconstructed jet can successfully be corrected. This also is valid for 
all centrality classes with different amount of background, as it is shown in figure |2|^b) 
for peripheral events (centrality 50 — 80%). 

The width of dpT is about lOGeV/c in the most central events (0 — 10%) and 
decreases for more pheripheral events (50 — 80%), where it is about 1 — 2GeV/c only. 
However, the Gaussian fit does not describe the distribution very well in such peripheral 
events with low number of tracks. 

In addition to the random cones from all (background) jets the 6pT distribution 
from the random cones without the two leading jets is shown (figj2| open circles). Here, 
the tail on the right-hand side almost disappears. This indicates that the tail indeed 
comes from some jets in the heavy-ion events and are not part of the general background 
contribution. 

The measured background, of course, directly depends on the chosen low p^ cut off. 
In the studies presented so far we benefit from the good tracking capabilities of ALICE 
for tracks of very low down to 150MeV/c. However, in view of jet reconstruction 
a higher p^ cut can help to reduce the influence of the background. For the track-px 
cut of 0.15GeV/c, which we used so far, is in centrality class — 10% the average 
p = 136GeV/c, it decreases to p = 61GeV/c for p^ > l.OGeV/c and further to 
p = 13GeV/c for p^ > 2.0GeV/c. But the main advantage is the smaller background 
fluctuations which we can expect due to the reduced number of tracks. The effect for 
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Figure 2: Background fluctuations 6pT in Pb-Pb collisions at y/s = 2.76 TeV for central 
(0 - 10 %) (a) and peripheral events (50 - 80 %) (b) with track pt > 150 MeV/c. 

those track-pT cuts is shown in figure |3j For > 2.0GeV/c is a = 3.2 GeV, while the 
mean of the distribution basically stays around zero. That means that the background 
subtraction method also works for higher p^ cuts. Altough, we enforce the bias on a 
hard fragmentation, and make the situation worse for the reconstruction of quenched 
jets accordingly, if we raise the pt cut. That we have to consider. 
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Figure 3: Background fluctuations SpT in central collisions (0 — 10 %) for track 
Pt > 0.15 GeV/c, px > 1-0 GeV/c and px > 2.0 GeV/c. 
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7. Conclusion 

We presented the first measurement of tlie background fluctuations in lieavy-ion events 
for cfiarged tracks with a low px cut-off of 150MeV/c. Based on the different methods 
random cones and fast embedding we see a good agreement with a background fluctuation 
of (7 lOGeV/c in the most central events. The subtraction of the average background 
is well under control also for events with large number of tracks. To reduce the 
fluctuations for the reconstruction of jets in a lower px region it might be necessary 
to introduce a higher track-px cut. 

8. Outlook 

We did not discuss our studies of the background fluctuations as function of multiplicity. 
The dependence of the number of tracks allows a comparison with the Poissonian limit. 
The background fluctuations are expected to approach with the Poissonian limit in case 
they are caused by uncorrelated sources only. 

But, also other effects contribute to a broadening of the fluctuations. E.g. we 
estimate the background fluctuations for different orientations to the event plane for 
this purpose. 

Furthermore we investigate the background fluctuations with embedding of full 
PYTHIA jet events and quenched jets (QPYTHIA and PYQUEN). With those probes 
our sample becomes more realistic, but also more biased according to the jet finder, so 
we observe jet splitting and have signs of back reaction. 
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